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ABSTRACT 

We test the theoretical prediction that the straightest dust lanes in bars are found in strongly barred 
galaxies, or more specifically, that the degree of curvature of the dust lanes is inversely proportional 
to the strength of the bar. The test used archival images of barred galaxies for which a reliable 
non-axisymmetric torque parameter (Qb) and the radius at which Qb has been measured (?"(Qb)) 
have been published in the literature. Our results confirm the theoretical prediction but show a large 
spread that cannot be accounted for by measurement errors. We simulate 238 galaxies with different 
bar and bulge parameters in order to investigate the origin of the spread in the dust lane curvature 
versus Qb relation. From these simulations, we conclude that the spread is greatly reduced when 
describing the bar strength as a linear combination of the bar parameters Qb and the quotient of the 
major and minor axis of the bar, a/b. Thus we conclude that the dust lane curvature is predominantly 
determined by the parameters of the bar. 

Subject headings: galaxies: kinematics and dynamics - galaxies: spiral 



1. INTRODUCTION 

Dust lanes in galactic bars have been observed in 
nearby galaxies for a long time. Photographic surveys 
from the beginning of 20 th century already allowed ob- 
servation of these features. Pease (1917) describes the 
prominent and nearly straight dust lanes in NGC 5383 
as a 'dark streak'. The same paper presents images of 
NGC 1068 in which a dust lane is clearly visible. The 
dusty nature of these features was identified at a much 
later stage. Sandage (1961) writes in his atlas that 'one of 
the major characteristics of the SBb(s) [galaxies] is the 
presence of two dust lanes leaving the nucleus, one on 
each side of the bar and extending into the spiral arms'. 

Dust lanes have been recognised as related to shocks in 
the gas flow in barred galaxies since Prendergast (1962). 
Such shocks are usually found in the leading edge of the 
bar, roughly parallel to its major axis. The location of 
these shocks also corresponds to the location of areas of 
high shear, which prevents star formation (Athanassoula 
1992; see Zurita et al. 2004 for a graphic illustration in 
the case of NGC 1530). Simulations show that the dust 
lanes are nearly straight and near the center of the bar 
when the galaxy has no inner Lindblad resonance (ILR) , 
which is known to be quite a rare occurrence. ILRs cause 
the dust lanes to be offset from the bar major axis, and 
to be curved with the concavity pointing to the bar ma- 
jor axis (Athanassoula 1992). Simulations by Patsis & 
Athanassoula (2000) show that the higher the gas sound 
speed, the smaller is the offset between the dust lane and 
the major axis of the bar. No dust lanes are expected in 
nuclear bars (Shlosman & Heller 2002). 

Athanassoula (1992) predicted from simulations that 
dust lanes would have a greater curvature in weaker bars 
and that dust lanes would be nearly straight in strong 
bars. This effect has been empirically verified using a 
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set of images of nine barred spiral galaxies by Knapen 
et al. (2002). The aim of the present Letter is to im- 
prove the statistics compared to the latter study, and to 
test definitively the prediction from Athanassoula (1992). 
Because dust lanes can be observed so easily they offer 
a fundamental handle on the underlying dynamics of the 
galaxy and the bar, yet they have hardly been studied 
observationally. A basic further aim of the present work 
is thus to explore in a statistical fashion the dependence 
of dust lane shape on the dynamical properties of the 
galaxy which hosts them. 

2. SAMPLE SELECTION AND BAR STRENGTH DATA 

We have selected all galaxies with a published non- 
axisymmetric maximum torque parameter for the bar, 
Qb, and a published radius at which Qb has been mea- 
sured, r(Qb). Its definition and the values for 266 
different galaxies can be found in Laurikainen & Salo 
(2002); Block et al. (2004), Laurikainen el al. (2004), 
Laurikainen et al. (2006), and Comeron et al. (2009). 
Some of these authors publish Q g , which includes the 
non-axisymmetries not related to the bars (but to spiral 
arms). However, Q g and Qb take very similar values in 
barred galaxies and can be interchanged for our purposes. 
In several cases Qb has been reported in more than one 
paper, in which cases we have used the most recent mea- 
surement (except for those in Comeron et al. 2009 which 
are based on shallower images). For NGC 1068 we pre- 
ferred the value from Comeron et al. (2009) to that from 
Laurikainen et al. (2004) because in the latter the au- 
thors measure Qb of the nuclear bar, whereas the dust 
lanes are hosted by the large bar. 

There are small variations in the techniques used to 
determine Qb in the various articles quoted above but 
we have checked that the Qb determinations are consis- 
tent within the error bars. These technical differences 
cause some scatter in the Qb measurements. The most 
important source of scatter is the use of different vertical 
scale-heights when measuring Qb and the large scatter in 
the observed vertical scale-height within one single Hub- 
ble type. For a detailed discussion see Laurikainen & 
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Salo (2002), where the authors consider that the uncer- 
tainties in both the orientation parameters and in the 
disc scale-height determination may cause an error of up 
to 10-15% in Qb- As we are mixing data from several 
sources we take a slightly higher estimate for the overall 
uncertainty in our Qb values of 20%. 

The 266 galaxies with a Qb determination in the lit- 
erature have morphological types spanning from SO/a to 
Sd. For those galaxies in the Sloan Digital Sky Survey 
(SDSS) DR7 (Abazajian et al. 2009) we examined single- 
band (g and i) and color-index (g — i) images to select 
galaxies presenting clear dust lanes. For all other galaxies 
we examined all optical images available in the NASA- 
IPAC Extragalactic Database (NED). Out of the original 
266 galaxies 55 have been found to host dust lanes. This 
is a lower limit to the true proportion of recognizable 
bar dust lanes in disc galaxies because for several galax- 
ies only Digital Sky Survey (DSS) images are available, 
on which bars are often saturated. The final sample for 
our dust lane study is thus of 55 galaxies. Table [1] lists 
these galaxies, their Qb, and its source. 

3. MEASUREMENT OF THE DUST LANE CURVATURE 

The curvature of the dust lanes was determined di- 
rectly from the SDSS and NED images using the proce- 
dure described by Knapen et al. (2002). We first depro- 
jected the galaxy images using the deprojection param- 
eters listed in Table Q] We then measured the change in 
the angle of the tangent to the dust lane in the range 
in which its curvature is constant, thus often ignoring 
the inner parts of the dust lanes where xi orbits in the 
circumnuclear region influence their shape and enhance 
the curvature (Athanassoula 1992). We also often ig- 
nored the outermost parts of the dust lanes, where they 
are closer to the corotation radius. The measure of dust 
lane curvature, Aa, as used here is arrived at by multi- 
plying the change in the tangent angle defined above by 
the only measure of the bar length that can be reliably 
derived in both observed and modelled bars, namely the 
radius at which the torque is maximal (r(Qb)). Aa is 
thus a dimensionless quantity which takes into account 
the size of the host bar. The results are listed in Table [TJ 

The uncertainty in the dust lane curvature is hard to 
quantify. The error due to galaxy deprojection is most 
probably not very large (< 10°) because in most cases 
dust lanes are intrinsically rather straight, and thus less 
sensitive to deprojection than a more strongly curved fea- 
ture might be. The uncertainty in determining the tan- 
gent to the dust lane is larger because the dust lanes are 
often not too well defined, and sometimes rather fluffy. 
By repeating our measurements for the same dust lanes, 
but using slightly different tangent points, we found that 
the results agree to within a few degrees per kpc. We 
thus adopt a conservative estimate of uncertainty in the 
dust lane curvature of <5Aa ~ 15°. 

4. RESULTS 

The relation between the bar strength as described by 
Qb, and the dust lane curvature, Aa, is shown in Fig.[TJ 
The plot has the axes flipped with respect to those used 
in Athanassoula (2002) and Knapen et al. (2002). The 
most striking feature in the plot is the clear division be- 
tween a populated lower left and an unpopulated upper 
right part, separated by a straight and well-defined enve- 




Qb 

Fig. 1. — Dust lane curvature related to Q^. 



lope. This implies that strong bars cannot have curved 
dust lanes, and that galaxies with smaller bar strengths 
can have progressively more curved dust lanes. A small 
unpopulated area appears in the extreme lower left, at 
the location of galaxies with straight dust lanes and small 
Qb- If on the vertical axis of the plot we use the change 
in the angle of the tangent to the dust lane (in °/kpc) 
instead of the dust lane curvature the results remain es- 
sentially the same. 

A similar diagram was presented by Knapen et 
al. (2002), but they only plotted nine data points. In 
their figure, the points appear roughly aligned, but in our 
new version, and thanks to the much larger number of 
data points, we can see how the relationship has a spread 
much larger than can be expected from the error bars. 
It is thus reasonable to suppose that although Qb limits 
the allowed dust lane curvature, it does not prescribe it. 
Other factors must determine the degree of curvature, 
especially for the weaker bars, and we used numerical 
simulations to explore which factors might contribute. 

5. SIMULATING DUST LANES IN BARRED GALAXIES 

We ran a set of 238 3D numerical simulations and stud- 
ied the dynamical response of the gas within a barred po- 
tential using the code FTM 4.4 (updated version) from 
Heller & Shlosman (1994). In each individual simulation, 
we used 10 5 isothermal, non-self-gravitating, collisional 
gas particles. We reproduced the potential of barred 
galaxies by means of a Miyamoto-Nagai disk, bulge and 
halo. The bar is reproduced by a Ferrers potential (Fer- 
rers 1877), which is introduced gradually, in one rotation, 
in a way that conserves the mass. 

We use different parameters for the bar ellipticity, 
mass, and pattern speed, and for the bulge mass ratio 
(B/B+D) in order to cover the complete range of galax- 
ies in our set of observations presenting dust lanes. 

The Miyamoto-Nagai potential is described by the 
equation: 



GM 



sJr 2 + (a + Vz 2 + b 2 f 



For the disk we use Ma 



5 x 10 10 M, 



0> 



(1) 



5.0 kpc 



and b — 0.5 kpc. For the bulge we use b = 0.05 varying 
B/B+D from 0.1 to 0.5 in steps of 0.1 by increasing the 
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TABLE 1 

Properties of the host galaxies and dust lanes in our sample. 
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Notes: galaxy ID (col. 1), distance to the galaxy in Mpc from Comeron et al. 2009 and LEDA (col. 
cllipticity from the literature (col. 4), dust lane curvature as defined in Sect. 3 (col. 5), Qj, (col. 6), 



for the galaxy disc PA, ellipticity, Q^, and r(Q^,) (col. 
al. (2004), and C09: Comeron et al. (2009). 

bulge mass. For the halo we use Mh — 10 11 
0.0 kpc, and 6—10 kpc. 
The Ferrers potential is described by 



2), galaxy disc PA (col. 3) and disc 
r (Qb) (col. 7), and literature source 



L02: Laurikainen & Salo (2002), B104: Block et al. (2004), L04: Laurikainen et 
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where m = - 
minor, and c the vertical axis, taken here in such a way 
that a/b takes the values 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5 ac- 
counting for the ellipticity of the bar. We use a — 6 kpc, 
c = 6/4 and two different values for the bar mass: 
M b = 5.3 x 10 10 AL Q and M b = 8 x 10 10 M Q . 

The pattern speed is chosen so that in a linear approx- 
imation the corotation radius is always greater than the 
length of the bar. In our set of simulations this means 
10, 20, 30, and 40 km sec" 1 kpc" 1 . 

As dust lanes are found in areas with gas shocks (Pren- 
dergast 1962), we inferred dust lane curvature in the sim- 
ulated galaxies by measuring the curvature of gas density 
enhancements using the same methodology as used for 
the dust lanes in real galaxies. This measurement was 
applied after two bar rotations. Well defined measure- 
able dust lanes appeared in 88 out of 238 simulations. 
The results can be seen in the left panel of Fig. [21 where 
the symbol coding indicates the ratio between the coro- 
tation radius, tcr, and the bar radius, tv For conve- 



nience, we define in the context of this Letter that for 
r CR < 3.0 rb we have a 'normal' bar and for 3.0 r b < rcn 
we have a 'very slow' bar (more details on this bar clas- 
sification can be found in Debattista & Sellwood 2000). 
rcn was always measured using a linear approximation 
and n, = 6 kpc was fixed in the simulations. 'Very slow' 
bars (two out of 88 in the simulations which present dust 
lanes) are rare in nature (just one is found in the compila- 
tion by Rautiainen et al. 2008) and thus are not included 
in the discussion in the next section because they have a 
non-standard behaviour. 

6. DISCUSSION 

The observational part of our study confirms that the 
degree of curvature of dust lanes decreases with increas- 
ing bar strength, i.e. stronger bars have straighter dust 
lanes (Fig. [T]). This was predicted by Athanassoula 
(1992) and confirmed in a preliminary way from observa- 
tions of nine galaxies by Knapen et al. (2002). Here, we 
find that the relation, in fact, outlines an upper limit to 
the dust lane curvature allowed for a particular value of 
Qb, in the sense that a bar with low Qb can have cither 
straight or curved dust lanes, but bars with high Qb can 
only have straight dust lanes. 

We performed a series of numerical simulations to ex- 
plore which parameters may be at the origin of the spread 
of the data points in the left panel of Fig. [2] We tested 
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0.0 0.2 0.4 0.6 0.8 -0.1 0.0 0.1 0.2 0.3 0.4 

O b Q b -0.087a/b 



Fig. 2. — Left panel: dust lane curvature related to Q b in the simulated galaxies. Filled dots stand for 'normal' bars (rcR < 3.0 r b ) 
and empty dots stand for 'very slow' bars (r^R > 3.0rt,). Right panel: dust lane curvature related to a combination of Qj, and a/b which 
minimizes the spread. The line is the best fit to the 'normal' bar data. There is no error bar on Q because it has been analiticaly calculated 
from the potentials used in simulations. 



whether this spread is caused by different bulge masses, 
different galaxy masses inside the corotation, and differ- 
ent bar pattern speeds. We find no effect of these pa- 
rameters on the spread except for the bar pattern speed: 
the spread for slow bars is higher but this effect is quite 
small. We found, however, that the data can be orga- 
nized in a plot with much smaller spread by fitting the 
logarithm of the dust lane curvature to a linear combi- 
nation of Qb and a/b which are parameters that depend 
only on bar properties. The best fit is 

Q b -0.087a/6= (0.156±0.020) - (0.119±0.015) log 10 Aa 

(3) 

and it is represented in the right panel of Fig.O The cor- 
relation coefficient of the fit is p = 0.66, better than when 
just fitting Qb and log 10 Aa (p = 0.58). This fit shows 
that Qb overestimates the effect of the bar ellipticity on 
the dust lane curvature. 

The correlation between the linear combination of Qb 
and a/b and the logarithm of Aa cannot yet be easily 
tested with observational data. Measuring the elliptic- 
ity of the bar requires an accurate decomposition of high 
quality images of the galaxy into its different compo- 
nents. Further work is needed in order to obtain these 
decompositions and refine the comparison with the sim- 
ulations. 

7. CONCLUSIONS 

Using 55 galaxies with published measurements of the 
bar strength indicator Qb, and for which we could obtain 
accurate measurements of the dust lane curvature, Aa, 



we confirm the theoretical prediction made by Athanas- 
soula (1992) that the dust lane curvature anticorrelates 
with the bar strength. Strong bars thus host straight 
dust lanes, while weaker bars can host more curved dust 
lanes. We do find that the anticorrelation has a large 
spread, but from a set of 238 numerical simulations of 
barred galaxies we show that this spread can be greatly 
reduced by using an appropriate linear combination of 
bar parameters (Qb and the axis ratio a/b). 
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